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Abstract. We consider the equations governing incompressible, viscous fluids in three space dimen- 
sions, rotating around an inhomogeneous vector B(x): this is a generalization of the usual rotating 
fluid model (where B is constant). We prove the weak convergence of Leray-type solutions towards 
a vector field which satisfies the usual 2D Navier-Stokes equation in the regions of space where B is 
constant, with Dirichlet boundary conditions, and a heat-type equation elsewhere. The method of 
proof uses weak compactness arguments. 



Resultats de convergence faible 
pour des equations des fluides tournants non homogenes 

Resume. On considere les equations modelisant des fluides incompressibles et visqueux en trois 
dimensions d'espace, en rotation rapidc autour d'un vecteur non homogene B{x): on generalise ainsi 
le modele habituel des fluides tournants (oii B est constant). On montre la convergence des solutions 
de Leray vers un champ de vecteurs qui verifie les equations habitucUes de Navier-Stokes 2D dans les 
regions de I'espace oii B est constant, avec des conditions aux limites de Dirichlet, et une equation de 
type chaleur ailleurs. La methode de demonstration repose sur des arguments de compacite faible. 

1. Introduction 

The aim of this article is to study the asymptotics of solutions of rotating fluid equations, in 
the case when the rotation vector is non homogeneous. We consider a domain = Jl/i x O3, 
where O/j denotes either the whole space R^ or any periodic domain of R^ , and similarly Q.^ 
denotes R or T. We are interested in the following system: 

dtu + u-Vu — uAu + -u A B + Vp = on R"*" x fi, 

e 

(1-1) V- n = on R+ x n, 

U|j^Q = on J7 

where B = be^ is the rotation vector, and 6 is a smooth vector field defined in 0,h- We shall 
suppose throughout this paper that b does not vanish, and is equal to a positive constant 60 
perturbated by a smooth function; more assumptions on b will be made as we go along. 
Before stating the result we shall prove here, let us recall some well-known facts in the constant 
case {b = 1). The rotating fluid equations, with b constant and homogeneous, modelize the 
movement of the atmosphere or the oceans at mid-latitudes (see for instance (HI or ^El)- The 
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fluid is supposed to be incompressible (which corresponds to the hydrostatic approximation), 
and its viscosity is > 0. The vector field u is the velocity and the scalar p is the pressure, 
both are unknown. The parameter e is the Rossby number, and its inverse stands for the 
speed of rotation of the Earth. Taking the limit e — > means that the scale of motion of the 
fluid is much smaller than that of the Earth. Note that one can also see S as a magnetic 
fleld, in which case it makes sense to understand what happens when b is not homogeneous; 
that also holds if one wants to study the movement of the atmosphere in other regions than 
mid-latitudes. 

In the constant case, those equations have been studied by a number of authors. We refer 
for instance to the works of A. Babin, A. Mahalov and B. Nicolaenko [21- I- Gallagher [Jj, 
E. Grenier JOl for the periodic case, and J.-Y. Chemin, B. Desjardins, I. Gallagher and E. 
Grenier for the whole space case as well as [H] for the case of horizontal plates with Dirichlet 
boundary conditions (for such boundary conditions we refer also to the work of E. Grenier 
and N. Masmoudi 11 as well as N. Masmoudi J^J- The results in those papers concern both 
weak and strong solutions; in this article we shall only be concerned with Leray-type weak 
solutions ([l2i): we will see in Section [2l below that their existence is an easy adaptation of the 
proof of Leray's existence theorem ^2j. In the constant case, it is known that weak solutions 
converge towards the solution of the two-dimensional Navier-Stokes equations. Such a result 
in the whole space case is due to Strichartz-type estimates (which are obtained by writing the 
solution of the linearized problem in Fourier space), whereas in the periodic case it follows 
from the study of the (discrete) spectrum of the rotating fluid operator (following methods 
introduced by S. Schochet in J7]). 

The problem here if we want to follow those methods is that it does not seem a good idea to 
take the Fourier transform of the operator 

Lu'^= P{uAB), V-u = 0, 

where P denotes the Leray projection onto divergence free vector fields, when B is not ho- 
mogeneous; moreover the study of the spectrum of L is not an easy matter. So our strategy 
to study this problem is first to try and recover the well-known results of the constant case 
without using any information on the spectrum of L (other than the determination of its 
kernel), and without using the Fourier transform. This will be achieved in Sectional Then 
the study of the variable case will be an adaptation of the constant case, in Section^ 

Before stating the results we shall prove in this paper, let us comment on the difficulties 
compared with the constant case: as stated above, it is easy to construct a bounded family 
of weak solutions to our problem, whether b is constant or not. Hence one can construct a 
weak limit point u, and the question we want to address is to find the equation satisfied by u. 
Of course the problem consists in taking the limit in the non linear part of the equation. As 
noted above, we do not wish to study the spectrum of the operator L since that seems to be 
a difficult issue. So we cannot apply the usual, constant b methods, as to our knowledge they 
all involve spectral properties of L. The idea therefore is to turn to what is known as "weak 
compactness methods", in the spirit of Lions and Masmoudi |13j-|14j (for the incompressible 
limit). We shall recall briefiy below what those methods are, and then we shall state the main 
results of this paper. 
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1.1. Weak compactness methods. Let us explain what weak compactness methods are all 
about. The idea is as follows: as usual the trouble to find the limit of the equation comes 
from the bilinear terms. They can be separated into three categories: 

• products involving only elements of the kernel of the penalization L, which can be 
shown to be compact (see Corollarv 12.5(1 : 

• products of elements of the kernel against elements of (KerL)^, for which one can take 
the limit since elements of (KerL)^ converge weakly to zero (see Corollarv 12. 5() : 

• products involving only elements of (KerL)-*-, which are the problem. 

The idea now is to prove that in the last situation, the limit is in fact zero for algebraic 
reasons: in previous works on rotating fiuids, that result was proved essentially by writing the 
product of two elements of (KerL)-*" by projection onto eigenvectors of L. In the periodic case, 
a "miracle" in the formulation yielded the result (see [2]- [3] or whereas in the whole space 
case, Strichartz estimates did the job (and the convergence was strong), see 0. In this paper 
we will show that the result has in fact not much to do with spectral properties of L, but is 
due to simple algebraic properties. Let us recall the result in the case of the incompressible 
limit, where such properties were first used (see jl3j). 

Proposition 1.1. [TSl Let (pe), (ue), (Oe) be bounded families of L'^{[0,T], H^{n)) such that 

Pe ^ p, Ue ^ u, 6e ^ 9 as e ^ 0. 

Assume that 

edtPe + V-Ue = 0, 
edtUe + V{pe + Be) = Se, 

2 

edtOe + -V- Ue = s^, 
where Se, 4^0 in L^{[0, T], iJ-'*(17)) for some s > 0. Then 

PV- {ue (S) Ue) PV- (n (g) u) and V- (ueOe) V- {u9) 
in the sense of distributions. 



Proof. This result has to be compared with the so-called "compensated compactness" the- 
orems, in the sense that the convergences of some quadratic quantities in pe,Ue,6e are es- 
tablished under the assumption that some combinations of the derivatives of these functions 
converge strongly in time to 0. The proof consists in checking that the acoustic oscillations do 
not bring any contribution to the limiting terms. We introduce the following decompositions: 



so that 



Ue = PUe + V'll^e, and Oe = ^^^—-^ +TTe, 



Pue and — ^ are bounded in VF^'i([0,r],i/-" (17)), and 



edtV^Pe + Vvr, = {Id - P)se ^ 0, eStvr, + ^AiPe = ^4^0 in ^^([0, T],H-'{n)). 

We shall note in the following Se {Id — P)se and S'^ —s'^. The incompressibility and 
Boussinesq relations 

V- n = 0, V(p + 9) = 
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allow to identify the limits 

Pue u, ~ -^9 in l2([0, T] X n), 

5 

V'i/'e^O, TTs ^0 inw - L'^{[0,T] xn), 
from which we deduce that, in the sense of distributions 

Ue®Ue — Vlpe (8) Vipe ^ U U, 
6eUe - TTeWe ^ Ou. 

The key argument is therefore the following formal computation (which can be made rigorous 
by introducing regularizations with respect to the space variable x) 

= lP{-dt{e7r,VA) + TTeSe + S'eVA), 

which shows that the contribution of the acoustic oscillations is negligible. □ 

Inspired by the previous computation, we shall in this article try to use a similar method 
in the case of rotating fluids: we refer to the proofs of Propositions 13.41 and 14.41 for precise 
computations. 

1.2. Main results. Since we consider incompressible flows, we introduce the following sub- 
spaces of L^(r2) and H^{i}) 

R = {ue L^{n)/V- u = 0}, Y = {u£ H\Q)/V- u = 0}. 

We will also use the following notation for the inhomogeneous Sobolev spaces 

H'{^) = {n G V'{^) I {Id - A)"/\ G L^iSl)}. 

Similarly homogeneous Sobolev spaces will be defined by 

R\Vl) = {n G V\Vl) I (-A)^/2^ G L^iSl)}. 

It will appear clearly in the following that the horizontal variables play a special role in 
this problem. Consequently we shall use the following notation: if x is a point in 17, then 
we shall note its cartesian coordinates by (xi,X2,X3), and the horizontal part of x will be 

denoted x^ ^= (xi,X2) G Similarly we will denote the horizontal part of any vector 

field f hy fh, the horizontal gradient by V/^ {di,d2) and its orthogonal by = {82, —di), 

and the horizontal divergence and Laplacian respectively by div/j/ difi +82/2 = V/^ • fh 

and Ah=^df + dl 
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Finally as usual, C will denote a constant which can change from line to line, and Vp will 
denote the gradient of a function which can also change from line to line. 

Before stating the main theorems of this paper, let us give some additional definitions. We 
will note by 

5 =^{x G 0/V&(x) =0} and O =^ {x G Q / ^ 0}. 

Finally S will be the interior of the singular set S. We will assume in the sequel that 
{HQ) The set Vl\{S U O) is of Lebesgue measure ; 

{HI) 5 is a smooth domain. 

On each connected component Oj of O, there is a smooth function cTj 
{H2) such that {b,aj,X3) is a global smooth coordinate system 

and Oj = {xh G | ib{xh),aj{xh)) G Bj x Sj)} x n^. 

Now we are ready to state the main theorems of this paper. The first result, rather standard, 
shows that there are weak solutions to the system 

Theorem 1. Let be any vector field in H. Then for all e > 0, Equation has at 

least one weak solution Ue G L°°(R+, H) n L^(R+, i^T^). Moreover, for all t > 0, the following 
energy estimate holds: 

(1.2) \\ue{t)\\l2+2iy [ \\Vue{s)\\l2ds <\\u°\\l2. 

Jo 

Now the aim of the paper is to describe the limit of ti^ as e goes to zero. We will first 
concentrate on the constant case. 

Theorem 2. Suppose that B = be^ where b is constant and homogeneous. Let vP he any 
vector field in H, and let Us be any weak solution of in the sense of Theorem CI Then 
converges weakly in L^^^(R^ x Q,) to a limit u which if = R is zero, and if Q3 = T is the 
solution of the two dimensional Navier-Stokes equations 

{NS2D) dtu - i^AhU + Uh ■ VhU = {-VhP, 0), divhUh = 0, U|t=o = / u°{xh,X3) dx^. 

Jt 

Remark 1.2. This theorem is by no means a novelty, it is even rather less precise than other 
such results one can find in the literature fjSj-jl], 0,0; CO])- we will see in Section\^ the 

interest of this result lies in its proof, which contrary to the references above, does not depend 
on the boundary conditions (which can be the whole space or periodic, in each direction). 

Now let us state the new result of this paper, concerning the case when b is not homogeneous. 

Theorem 3. Suppose that B = be^ where b = b{xh) is a nonnegative smooth function, say 
a C^(J7ft) perturbation of a constant, and where assumptions {HO) to {H2) are satisfied. Let 
be any vector field in H, and let be any weak solution of i|I._ZI) in the sense of TiieoremQJ 
Then converges weakly in Lf^^{'R.~^ xQ) to a limit u which if = R is zero, and if O3 = T 
is defined as follows: the third component satisfies the transport equation 

dtlf - i^AhV^ + Uh ■ VhU^ = 0, ^3^1^ = 0, u^^Q = j u^'^{xh, X3) dx3 in R+ x (O U S), 

while the horizontal component depends on the region of space considered: 
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• in S, Uh satisfies the two dimensional Navier-Stokes equations {NS2D) with Dirichlet 
boundary conditions. 

• in O, Ufi satisfies the following heat equation: 

dtUh - i^UAhUh = 0, 

where H is the orthogonal projection onto the kernel of L (which can be extended 
to T>'[0 U 5) J which satisfies in particular 

- {IlAhUh\uh)L2(^0) = \\^hUh\\\2f^Qy 

Remark 1.3. In the regions where b is homogeneous, we recover at the limit the 2D Navier- 
Stokes equations as usual. The Dirichlet boundary conditions appear quite naturally, consid- 
ering that on the other side of the boundary one finds that is proportional to V^6 which 
vanishes on the boundary of S. More surprising is certainly the result in the region where b 
is not homogeneous. This can be understood as some sort of turbulent behaviour, where all 
scales are mixed due to the variation of b. Technically the result is due to the fact that the 
kernel of L is very small as soon as b is not a constant, which induces a lot of rigidity in the 
limit equation. 

The structure of the paper is as follows. In the next section, we present the operator L and 
study its main properties (proof of Theorem ^ kernel, projections onto KerL). The following 
section is devoted to the proof of Theorem [2 Although the result is not new, we present an 
alternative proof which holds regardless of the domain (with no boundary). This serves as a 
warm-up to the final section, in which the general variable case is presented, with the proof 
of Theorem 131 

Remark 1.4. One can wonder about what remains of those results under more general 
assumptions on the rotation vector B. First let us consider the case when B = b{t,Xh)e^ 
depends also on time. The results of Section are identical, and if the sets S and O are 
independent of time one recovers the same type of theorem as in the constant case. In 
particular the equation on u in O is derived in an identical way to Section \4.2l If the sets S 
and O do depend on time, then one has to be a little bit more careful and this issue will not 
be treated here. A more physical problem is the case when the direction of B is not fixed, 
in other words when B is a three component vector, depending on all three variables. Then 
geometrical problems appear, simply to determine the kernel of L ; this will be dealt with in 
a forecoming paper. 

2. Study of the singular perturbation 
2.1. Energy estimate. In this section we shall prove Theorem ^ stated in the introduction. 

Proof. The structure of the equation (|1.H) governing the rotating fluids is very similar to 
the one of the usual Navier-Stokes equation, since the singular perturbation is just a linear 
skew-symmetric operator. Therefore weak solutions " a la Leray " can be constructed by 
the approximation scheme of Priedrichs : approximate solutions are obtained by a standard 
truncation J„ of high frequencies. In order to obtain uniform bounds on these approximate 
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solutions, we have just to check that the energy inequahty is stiU satisfied. Computing formally 
the scalar product of by u leads to 

1 d „o ,0 . 1 

Integrating by parts (without boundary) and using the incompressibility constraint, we get 

which holds for any smooth solution of 



|u|||2 = — J(^2^u-\/)\u\'^ — vu-/\u H — u-u A B + u-VpJ dx. 



W\\l2 = -Z^||VU||^2, 



The energy inequality for weak solutions is obtained by taking limits in the approximation 
scheme. □ 



In particular, the energy estimate provides uniform bounds in L°^(R+,H) n L^(R+,iJ^) on 
any family {us)e>o of weak solutions of (jl.ll) provided that the initial data belongs to H. 

Corollary 2.1. Let be any vector field in H. For all e > 0, denote by a weak solution 
of i^nj)- Then there exists u E L°°(R~'",H) n L^(R^,i7^), such that, up to extraction of a 
subsequence, 

Uc —^u in w-Lfg^{'R'^ x ft) as e —)■ 0. 



2.2. Characterization of the kernel. We are interested in describing the asymptotic be- 
haviour of (tie), i-6- in characterizing its limit points. Of course, the equations satisfied by 
such a limit point u depend strongly on the structure of the singular perturbation 

(2.1) L : -u G H P(n AS) G H 

where P denotes the Leray projection from L^{fl) onto its subspace H of divergence- free 
vector fields. In particular, we will prove that u belongs to the kernel Ker(L) of L, which is 
characterized in the following proposition. 

Proposition 2.2. DeEne the linear operator L by |2Jj). Then u G H belongs to Ker(L) if 
and only if there exist V^c^ G L?'{fth) and a G L?'{Q.h) with 

such that 



Proof. We have 

P{uAB) = 0. 

Then, in the sense of distributions, 

rot {uAB)= 0, 

which can be rewritten 

(V- B)u + {B-V)u- {u-V)B - (V- u)B = 0. 
As V- = V- ti = and B = 663 , we get 
(2.2) bd^u- {u-V)be-i = Q. 
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In particular, dsui = d^U2 = from which we deduce that 
(2.3) ui,U2(^L^{Vlh). 

Note that in the case where VL^ = R, the invariance with respect to X3 and the fact that u G 
L^(r2) imply that ui = M2 = (and therefore us = by the divergence free condition). 

Differentiating the incompressibility constraint with respect to X3 leads then to 

^33^3 = -dl^ui - dl^U2 = 
in the sense of distributions. The function d^u^ depends only on xi and X2, and satis- 
fies / d'iUsdxs = 0. So d'iU-i = and 



(2.4) U3 G L^Qh), diui + d2U2 = 0. 

Combining 1)2. 3|) and 1)2. 4j) provides the existence of VhV G L^i^h) such that 

Ui = d2(p, U2 = —Blip. 

Replacing in (|2.2I) leads to 

Vh^-Vhb = bdsus = 0, 

which concludes the proof. □ 

Before applying this result to the characterization of the weak limit u, let us just specify it in 
two important cases. If V& = almost everywhere, u G H belongs to Ker(L) if and only if 

u = V^((? + ae3, 

for some V^(/? G L'^{^h) and a G L'^{^h)- If V& 7^ almost everywhere (in other words, ifQ\0 
is of Lebesgue measure zero), then the condition arising on u is much more restrictive : u G H 
belongs to Ker(L) if and only if on each connected component of O, 

n = F(6)V^6 + ae3, 
for some square integrable function F{b) and some a G L'^{Qh)- 

From this characterization of Ker(L), we deduce some constraints on the weak limit u. 

Corollary 2.3. Let be any vector field in H. Denote by (ne)e>o a family of weak solu- 
tions of i)I.I|) . and by u any of its limit points. Then, there exist ip G L'^ {Yi'^ , H"^ {0.^)) n 
L2(R+,ij2([7^)) and a G L~(R+,L2(1^;,)) nL2(R+,ii'i(0/,)) with 

such that 

u = V^ip + aes- 

Proof. Let x ^ I'(R^ x Q) be any divergence-free test function. Multiplying ()l.lj) by ex and 
integrating with respect to all variables leads to 

Ue{edtX + E'Ue'Vx + £^^X + X ^ B)dxdt = 0. 
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Because of the bounds coming from the energy estimate, we can take hmits in the previous 
identity as e — > to get 

J u AB-xdxdt = 0. 

This means that there exists some p such that 

u A B = Vp. 

As Ug satisfies the incompressibihty relation for all e > 0, 

V- M = 0. 

Then u{t) G Ker(L) for almost all t G R"*", and we conclude by Proposition 12.21 that there 
exist (f G L°^{'R+,H'^{nh))nL'^(R+,H^{nh)) and a G L°°{'R+ ,L'^{nh)) D L'^(R+ ,H^{nh)) 
with 

such that 

□ 

2.3. Decomposition by projection on the kernel. To go further in the description of the 
asymptotic behaviour (i.e. in the characterization of u), we have to isolate the fast oscillations 
generated by the singular perturbation L, which produce "big" terms in but converge 

weakly to 0. 

Therefore we introduce the following decomposition 

Ue = Us+ We, 

where = nn^ is the L? orthogonal projection of on Ker(L) and Wg = n_|_n£ is the 
projection of Ug on Ker(L)-'-. 

We have seen in the previous paragraph that the characterization of the kernel Ker(L) is 
strongly linked to the geometry of the vector field b. In order to obtain further regularity 
properties on 11 and n^, we then need a precise description of the singular set 

S = {x en/Vb{x) = 0}, 

which justifies Assumptions {HO) to {H2) given in the introduction. 

Before stating the main properties of 11 and Il±, let us give the following definitions: by 
Assumption {HO) it is enough to describe the limiting function u on O L) S. So it is natural 
to define the following function spaces: for all s > 0, H''{0 U S) is the closure of C^{0 U S) 
for the H^ norm, and we will note, for s >0, H~^{0 U S) the dual space of H'^{0 U S). 

It will be useful in the following to note that 

(2.5) Vs>0, H'{0US) C H'{n) and Vs < 0, H'{n) C H'{0 U S). 

Proposition 2.4. Define the linear operator L by J|) . Denote by H the orthogonal projec- 
tion of H onto Ker(L) and by Il± = Id — Il the orthogonal projection of H onto Ker(L)^. The 
operators H and Il± so defined have natural extensions to tempered distributions on O L) S, 
and for all s G R, there exists some Cs > such that for any function u G H^{0 U 5), 

I|nn||/^s(ciu5) < Cs||n||j^a(ciu5) and ||n^u|l//s(ciu5) < Cs\\u\\Hsi^Qyjsy 
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Proof. By Proposition I2.2( for all u € H 

Uu = V^99 + aes, 
for some Vji^ G L?'{ilh) and a G Lp'{Q.h) with 

By definition n_|_it = u — Tiu is orthogonal to any element of Ker(L). So for all {3 € L^(J7/i), 

(n — nii)-/3e3(ix = 0, 



which implies that 

(2.6) a = 1^ j usdxs. 

In order to determine ip, we consider separately the domains O and S. 
On 5, (/? is defined as 

"^hf = ( 77^ / '^hdx; 
and the smoothness properties stated in Proposition 12 .41 are obvious. 

On O, we use Assumption {H2) which implies that on each connected component Oj of O, u 
can be written u{x) = fj{b,aj,X3) where fj has the same smoothness as u since the change 
of coordinate is in C°°{Q). Then clearly the projection 11 is simply defined by 



Hit =^ 1 [ [ fj{b,aj,X3)dajdxs, 



^j\ JEj Jq3 

and the result follows. □ 

Corollary 2.5. Let be any vector field in H. Denote by (Me)£>o a family of weak solutions of 
and hyu any of its limit points. Consider a subsequence of (us) (abusively denoted (ue) ) 
such that 

Ue in w-Lf^^{Y{,'^ x $7) as e ^ 0. 
Then, ifu^ = Hus is the projection of on Ker(L), 

Ue strongly in Lfo^(R+ x (O U 5)) as e ^ 0. 

Proof. By Proposition l2.41 the projection 11 is continuous in L^. Then, by the energy estimate, 

Me is uniformly bounded in Lf^jYi'^ ,}C) 
for some compact subset /C of which provides regularity with respect to space variables. 

The second step consists in getting regularity with respect to time. Apply 11 to the convection 
equation in : 

(2.7) dtUe + Il{Ue-VUe - v/\Ue) = 0. 

As Ue is divergence-free, 

Ue-VUe = V- (tie (8) U^) is uniformly bounded in L°°(R+, ^^"^'^(Sl)). 
From the energy estimate, we also deduce that 

Aiie is uniformly bounded in L^(R^, i7^"'^(r2)), 
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from which we infer that 

Ue-Vue — uAue is uniformly bomided in Lf^^lIV^ , H^^^^^ {Q)). 
Combining this with ()2.5p . we get by Proposition 12.41 

dtUg = — n(M£-Vu£ — i^Aug) is uniformly bounded in hj^^iji.^ , H^^J'^ {O U5)), 

which provides the expected regularity in t. 

Aubin's lemma ^ then gives the following interpolation result 

(us) is strongly compact in Lf^^iYi^' x [OVJ S)). 

By Proposition im H belongs to C {L"^ {^) , L'^ {il)) and therefore to C{w - L'^{n),w - L'^{^)), 
from which we deduce that 

Us = IlUe ^ I\E = U. 

Combining both results shows that Ue converges strongly to u in Lf^^{'R'^ x (O U 5)). □ 
2.4. Remarks concerning the regularity. 

2.4.1. Comparison with the gyrokinetic approximation. As mentionned in the introduction, 
the study of the asymptotics for an inhomogeneous penalization is a natural question in the 
magnetohydrodynamic framework, when B represents the magnetic field. Such a study has 
been performed for the gyrokinetic approximation 8^, that is for a kinetic model perturbed 
by a singular magnetic constraint : 

• in the case where B = 5(2:^)63, the singular limit is exactly the same as in the constant 
case : the fast rotation has an averaging effect in the plane orthogonal to the magnetic 
lines; 

• in the case where B has constant modulus but variable direction, extra drift terms are 
obtained due to the curvature of the field. 

A simplified version of this result can be written as follows. 

Theorem 4. |BJ Let be a function of L^{Q x R'^), and (/e) be a family of solutions of 

dtfe + v-VJe + ^v AB ■ V„/e = 0, t G R;, {x,v) G x R3, 
with initial condition 

fe{\t = 0) = f. 

Then the family {f^) is relatively compact in w * —L°°{^j^ x x R'^), as well as the family 
(gs) defined by 

9eit, X, w) = fe{t, X, R{X, --)w) 

where R{x,6) denotes the rotation of angle 6 around the oriented axis of direction B[x). 
Moreover, 

• if B = he^ with b G C^{Q.h,'i\,*^_), any limit point of {ge) satisfies 

dtg + vsd^^g = 0; 
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• if B G C^{Q) with Vx ■ B = and \B\ = 1, any of its limit points satisfies 
dtg + {w ■ B)B ■ V^g = \w A {3{w.B){B A VbB) -BA V^B - VbavjB) ■ V^g 

def 

with the notation Vy^ = ^ • V<I>. 

The result obtained in this paper is very different because of the incompressibihty constraint, 
which imposes a lot of rigidity to the system. In particular, the kernel of the penalization is 
much smaller and the limiting system has less degrees of freedom. 



2.4.2. A remark in the inviscid case. The weak compactness method used here allows to 
study the singular limit without regularity with respect to the time variable. However it uses 
crucially the strong compactness in x given by the energy estimate (|1.2|) . Implicitly we have 
actually considered the penalization 

Le:u£Y^P{uAB)- sAu € H-^{n). 

That rules out the possibility to manage an analogous study for inviscid rotating fluids, the 
first obstacle being to prove the existence of solutions for 

dtUe + {Ue-V)Ue + ^Ue A B + Vp = 0, V- Ue = 

on a uniform time interval [0, T]. Indeed the operator exp{tL/e) is not even uniformly bounded 
on H'{n) for s > ^■ 

Lemma 2.6. Define the hnear operator L by Then, the group (exp(tL))(gR generated 

by L is not uniformly bounded on H'^{Q,) for s > 

Proof. The proof of that result is simply due to the fact that by definition of 11 seen above, 
the trace of IIu is not defined on dS even if n G V. So 11 is not continuous on H^{^}) for 
s > i- Then the formula 



1 



n= lim - / e^^dt 



T 



implies that e*-^ cannot be uniformly bounded in H^{^) for s > i which proves the lemma. □ 



3. The case of a constant vector field B : 
THE 2D Navier-Stokes limit 



In the previous section, we have obtained a constraint equation on the limiting velocity field, 
which expresses that u belongs to the kernel of the singular perturbation L. This comes 
from the fact that the component We = 11^ has fast oscillations with respect to time, and 
consequently converges weakly to 0. In the case where = R, this characterizes completely 
the weak limit n = 0. 

Then it remains to get an evolution equation for u in the case where ^3 = T. A natural 
idea consists in projecting the evolution equation (jl.lj) for on the kernel of L, and to study 
its limit as e — s- 0. The difficulty is to take limits in the nonlinear terms : as Corollarv 12.51 
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provides strong compactness on the non-oscillating component , the problem comes actually 
to prove that the oscillating terms We do not product any constructive interference. 

In order to have a good understanding of this phenomena and of its mathematical formulation, 
we propose to consider first the case where the vector field B is constant and homogeneous. 
The convergence result established here is not so precise as the ones given in [2]-Il]) [Z| or jlUj . 
since it does not describe the oscillating component and consequently does not provide any 
strong convergence. Nevertheless the proof is less technical (in particular it does not require 
any knowledge on the spectral structure of L), which allows to consider more general cases in 
the sequel. 



3.1. Projection on the kernel. In order to obtain the evolution of the limiting velocity 
field n, the idea is to use the strong convergence 

Ue = Uue^u in lL(I^^ X (O U 5)), 

rather than the weak convergence 

Us —^u in w-Lf^^(R'^ x $7). 

Having this idea in view, we first determine the evolution equation for u^. 

Proposition 3.1. Let u^ be any vector field in H. For all e > 0, denote by u^ a weak solution 
of and by Us = Hu^ its projection on Ker(L). Then, 

(3.1) dtUe — n(ue A rot Ue) - v/Sms = 0. 

Proof. Identity ()3.1|) is essentially a variant of ()2.7|) . Indeed, in the case of a constant B, the 
projection 11 commutes with any partial derivative, in particular with the Laplacian A : 

n(z/AM£) = Z^AIIlie = Z^Alie. 

Then the key argument is the following identity : 

(3.2) u A rot n = V-^^ V- {u®u) + uV- u 

As is divergence-free, 

UeV- Ue = 0. 

As Ker(L) is embedded in the space of divergence-free vector fields, 

n(vM)^o. 

Replacing in (|2.7j) leads to the expected result. □ 

3.2. Brief description of the oscillations. 

Proposition 3.2. Let u^ be any vector field in H. For all e > 0, denote by u^ a weak solution 
of and by = Ii±Ue its projection on Ker(L)^. Then there exists We G Lf^^{'R'^ x Q) 
such that We = d^We- Moreover, 

edt {d^We - W,,3) + ^^We^B = re, 

^^■^^ edt{d2We,l - diWefl) + b{diWe,l + d2We,2) = Se, 

where Ve, Se converge to in Lf^jYl'^ , Hj^J'J'^ {O U S)) as e ^ 0. 
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Proof. As we have supposed in this section that S = Q, the projection of on Ker(L) 
satisfies, due to Proposition 12.41 



Then 



J Wsdxs = y ~ UUs)dx3 = 0, 



and there exists Ws such that 

We = dsWe- 

Moreover, as = T, we can always choose so that 

WedX3 = 0. 

It remains then to determine the equations for W^. Equation Hl.lj) imphes 

edtWs + W£ A be3 + Vp = euIl±Aus — eIl±{us-Vu£), 
which can be rewritten in terms of 

(3.4) edtdaW, + daW^ A he^ + Vp = S,, 

with 

Se = eiuIi^AUe - H^iUe-VUe)). 

From the energy bound, we deduce that the right hand side in 1)3.4(1 is of order e in the 

space -L^^^(R"'", Hj^^^'^{0 U S)). Indeed, using the continuity properties of n_L as weh as (|2.5|) 
we get 

l|n_LAne||i2([o,T],/f-i(C'U5)) < \\^'^e\\L-i([0,T],H-^(n)) < ll^e llL2([o,T],Hi(n)) 

and 

||n_L('Ue-VUe)||^2([o,T],/f-3/2(X)) ^ ll^e ''^^e llL2([o,T],H-3/2(i4')) 

< C||u£|jicx)(R+^i2(Q))||V'Ue||2,2([o,T],L2(n)) 

for all compact subsets K oi O yj S. Moreover, the right-hand side in (j3.4|) belongs to the 
image of Yl± and can therefore be written as a partial derivative with respect to X3, 

The equation on the third component provides then 

edtdj,We,'i + dsp = dsR^ 

where converges to in Lf^^(R^ , H^^^^'^ {O U S)) as e ^ 0. Integrating with respect to X3 
provides, since / p dx^ = 0, 

edtWe,3+P = Re- 

Replacing in (|3.4I) leads to 

edtd^We + d3We A 663 + V(i?, - edtWe,3) = Se 
which is the first identity in Proposition 13.21 
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In order to establish the second identity, we compute the rotational of t^'AA^ and write its last 
component 

Moreover, the right hand side in ()3.4|) belongs to the image of 11^ (recall that in this section n_|_ 
commutes with all derivatives) and can therefore be written as a partial derivative with respect 
to X3. Then integrating with respect to X3 leads to the expected equality, where the right-hand 
side converges to in Lf^^(R'^ , Hj^^^^"^ {O US)) as e goes to 0. □ 



3.3. Study of the coupling. The algebraic structure of the propagator (|3.3|) implies that 
the oscillating terms cannot interact and produce some contribution in the limiting equation 
governing u. Indeed the nonlinear term We A rot We can be rewritten as the sum of a total 
derivative with respect to X3 and a total derivative with respect to t/e, modulo a remainder 
which converges formally to 0. Then, in order to prove a rigorous convergence result, the first 
step is to introduce a regularization of the equations (|3.3|) and to get a control of the source 
terms in some strong norm. 

Lemma 3.3. Let be any vector field in H. For all e > 0, denote by a weak solution 
of inj) in L°°(R+, H) n L2(R+, ifi), and by = U±Ue its projection onto Ker(L)-L. Then, 
for all 6>0, there exists G Lf^jR+, n^i?^(r^)) such that 

(3.5) We - dsW^ in L;2^^(R+ x Q) as d ^ uniformly in e > 0, 

and 

edtid^W^ - VWl,) + dsW,' AB = , 

edt{d2wl, - diwl2) + Hdiwl, + 92^^,2) 

where for all 6, and sf converge to in Lf^^{lV~,Lf^^{^})) as e 




Proof. We introduce the following regularization: let k € C^(R^,R"'") such that k{x) = 
if \x\ > 1 and J ndx = 1, we define 



1 




and 

SO that = dsW^. 

By the energy estimate, for all T > 0, 

Ue is uniformly bounded in L'^{[0,T],H^{il,)). 
By Proposition 12.41 we infer that 

Ws = n_|_ne is uniformly bounded in L^([0, T],/C) 
for some compact /C of L^(r2). The result (|3.5|) then follows from the following fact: 
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hence there is some continuity modulus uj such that 

and the result follows. 
Regularizing (|3.,3j) leads to 

edt{d2Wl^ - d^wl^) + b{diwl, + d2Wl2) = s, * ks, 

because b is homogeneous. Then, for all T > and all compact subsets K oi OU S and for 6 
small enough, 

Ikf llL2([o,r],L2(K)) = Ike * I^5\\l-^{[0,T],L^{K)) 

< C'||K;5||^y5/2,l(R3)||re||^2([0,T],H-5/2(i^')) 

c 

where K' is a compact subset of O U 5 such that 

{x en/d{x,K) <d}c K'. 

And, in the same way, 

5 ^ 
lkellL2([0,T],L2(K)) < ^lkellL2([0,T],H-5/2(i^'))- 

For a fixed 5, Proposition gives the expected convergences. □ 

Equipped with this preliminary result, we are now able to study the coupling between the 
oscillating terms and to prove the following proposition. 

Proposition 3.4. Let he any vector field in H. For all e > 0, denote by a weak solution 
of ifl.Ip . and by = H_\_u^ its projection onto Ker(L)-'-. Then, 

pjweA rot Wedxz in D'(R+ x as e ^ 0, 

where P denotes the Leray projection. 

Proof. We start by proving that 
(3.7) 

P (^j We A rot Wsdxs — J /\ rot w^dx^ ^0 in X''(R^ x 0) as 5 — > uniformly in e. 
From the identity (|.S.2j) and the relations V- We = V-w^e = 0) 

we deduce that 

P (^We A rot We — We A rot We^ 



= -PV- {{We -W^^)(^We)- PV- {W° ® {We - <)) 

By Lemma EISl we get (|3.7j) . 
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It remains to prove that, for any fixed 6 > 0, 

(3.8) /\ rot w^^dx3 ^ as e ^ 0. 



Define = rot 




-d,Wl, + d3Wl, 

which is uniformly bounded with respect to e in L^{[0,T],nsH'{n)) by Lemma lOl We have 

\ dsWl, J \ dspl. 
From (|3.6() and the divergence- free relation d^W^^ = -diW^^ - SsW^^g' 

we deduce that the 

previous term can be rewritten 

eb^^dt{d3W^2-d2W^^ ^-^-^ 



e,3J " ' e,2 









d3pi2 




K dspU J 




or equivalently 

/ -eb-^du 

' ^^-'^ ■^e,2 

Integrating by parts with respect to X3 leads then to 

edtipl^dspU) + d^ipydtpU) + «i53pf,3 + s%pI^) 
1 ' edtipi^d^pl^) + d^ipydtpl^) + «253pf,3 - 4dwii) 
edt{pl,d3pl^) + d^ip^dtpy - «293pf,2 + rlidzpii) 
from which we deduce that, for any fixed (5 > 0, 

P A rot w^^dx3 



in the sense of distributions, as e ^ 0. 

Combining 1)3 .7^ and H3.8|) gives the expected convergence. □ 



3.4. Passage to the limit. In order to determine the limiting velocity field u, we have now 
to take limits in 1)3. If) which can be rewritten 

dfUir — Il{us A rot Us + A rot He + He A rot We + We A rot Ws) — i^Aus = 0, 

using the decomposition 

Ue = Ue+ We, 

where we recall that 

Us strongly in L^^^(R^ x 17), 
tt^e ^ weakly in L^„^(R+ x Vl). 
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Then standard arguments using (|H.2|) show that 

dtu — n(u A rot u) — v/SJI = Hm H{we A rot We), 

and by Proposition 13.41 we get 

dtu — n(n A rot u) — vHvd = 0. 

In the case where B is constant, the projection 11 reduces to a simple averaging with respect 
to X3, and we recover the usual convergence result towards the 2D1/2 Navier-Stokes equation. 

4. The case of a variable vector field B : 
a turbulent behaviour 

In this section we shall prove Theorem|31stated in the introduction, concerning the general case 
when the rotation vector B = h{xh)ez is inhomogeneous. We suppose that Assumptions {HQ) 
to {H2) are satisfied. If Q.^ = R, then u = simply because it is in L^(r2) but only depends 
on the horizontal variables. So from now on we can suppose that = T. 

The strategy of proof is quite similar to the constant case, so we shall often be referring to the 
results of the previous section. The first remark to be made is that if b is constant in some 
positive measured region of O/j, then in that region the results of the previous section should 
apply and one should recover at the limit the usual two-dimensional behaviour. 

Moreover, the results of Section [21 hold for any B, so in particular any weak limit point u 
of a sequence of weak solutions Ue to is in the kernel of L according to Corollary 12.31 
That means in particular that the third component does not see the difference between O 
and S since the elements of the kernel of L have the same third component whether h is 
homogeneous or not. So in the following, we shall restrict the study of the limit system to 
the horizontal components only. As in the previous section, the proof of Theorem |31 consists 
in finding the equation satisfied by u (at least its horizontal part u/j), by taking the limit of 
the equation satisfied by the horizontal part of Us = IIu^ . The first result we shall establish 
is that in the general, variable h case, there is no coupling between oscillating vector fields 
yielding extra terms in the averaged equation. This will be a generalization of Proposition E3] 
to the variable case, and the analysis will follow closely the proof of Proposition 13.41 Then we 
shall write the averaged equation on O. Finally we shall concentrate on the S case and show 
the limit Uh satisfies a two-dimensional Navier-Stokes equation with homogeneous boundary 
conditions on the boundary of S. 

4.1. The averaged equation. Let be a family of weak solutions to and define Ue = 

Hue- Recalling that the elements of Ker(L) are divergence free, we have as in the constant 
case 

dtUe — z/IIAue — n(iie A rot Ug) = 0. 

Of course the projector 11 no longer commutes with (horizontal) derivatives. However as 11 
belongs to C{w — H^{0 \J S),w — H^{0 U 5)) for s < 0, we clearly have as e goes to zero, 

IIAue ^ nAu. 

Let us now take the limit in the quadratic term. 
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Lemma 4.1. Let G L^(r2) be a divergence-free vector field. For all s > 0, denote by Ug a 
weak solution of and by Us = Hue (resp. Ws = n^Uej its projection onto Ker(L) (resp. 

Ker(L)-'-j. Then the following results hold in D'(R'*' x (O U S)), as e goes to zero: 

(4.1) U{ue ■ Vus) U{u ■ Vu), 

(4.2) U{w, ■ VUe + Ue ■ VWe) ^ 0, 



(4.3) Uiwe ■ Vwe) ^ 0. 

Proof. The results (|4.1() and 1)4. 2() are simply due to the compactness of Ue as well as the 
fact that Ws goes weakly to zero, results given in Corollarv 12.51 We also use the continuity 
properties of 11 stated in Proposition 12.41 The more difficult result to prove is of course 1)4. 3|) . 
The method will follow the proof of Proposition 13.41 and will be achieved in two steps. First 
we show that one can smoothen out the equation satisfied by Wg, and then we perform some 
algebra on the bilinear term in the equation, as in the constant case. We shall therefore 
continuously be referring to the methods of Sectional 

Let us start by proving the following result, analogous to Proposition 13.21 

Proposition 4.2. Let G L'^{Q) be a divergence-free vector field. For all e > 0, denote 
by a weak solution of and by = n_|_iie its projection onto Ker(L)-'-. Then there 

exists We^s G L^^^(R+ x 0,) such that Ws^s = 93VFe,3. Moreover, 

(4.4) 

edt{d2We,l - diWe,2) + diVh{bWe,h) = Se, 

where f^, converge to in Lf^^{R:^,H^^J'^{0 U S)) as e ^ 0. 

Proof. We shall omit the proof of that result here, as it is identical to the proof of Propo- 
sition 13.21 we just have to notice that the third component We^j, is of vertical mean zero, so 
can as in the constant case be replaced by d^We^'i- The other components, contrary to the 
constant case, cannot be transformed in that way, so remain as they are. The rest of the proof 
is identical to the constant case. □ 

Now as in the constant case, let us smoothen out Equation (|4.4|) . 

Lemma 4.3. Let G L^(J7) be a divergence-free vector field. For all e > 0, denote by Ue 
a weak solution of i^HJ); ^nd by Wg = n_|_tt£ its projection on Ker(L)-'-. Then, for all 6 > 0, 
there exists w^ f^ G Lf^^{Y{ + ,r\sH'{n)) and W/g G Lf^^{Yi+ ,r]sH'{Q)) such that 

(4.5) We,h — wlf^^ Q in Lfg^{K^ x 0) as 5 ^ uniformly in e > 0, 



(4.6) u'e^s - d'iWl^ in Lf„^(R+ x 9.) as 5 ^ uniformly in e > 0, 

(4.7) '^ll^/i^f,/illL2(n) — > as (5 — > uniformly in e > 
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and 

sdtiw' - VhW[,) + wt, AB = 



e 



(4.8) 

where and converge to in Lf^^{'R'^,Lf^^{0 U S)) as e and 5 go to 0. More precisely 
for any subset K of O U S there is a constant C (independent of e and 6) and a constant cs 
depending only on 5 such that 

(4.9) \\r', + 4h2iK)<C50e{l) + C5. 

Proof. We shall not write all the details of the proof here, since it is very similar to the 
constant case f Lemma 13. 3() : let us simply point out where the fact that b is not constant 
appears — note that in 1)4. 9() . the part QOe(l) is precisely due to the terms of the constant 
case, and we will see here that the fact that b is no longer constant yields terms which are 
estimated by C6. In the approximation of the equation, the only difference with the constant 
case is that of course {we A B) * ks is not equal to (wg * ks) /\ B. Moreover of course 1)4. 7() is 
obvious in the constant case, since n_|_ commutes with partial derivatives. So we need to deal 
with those two problems due to the fact that b is not constant. 

First of all, the difference between (ws /\ B) * ks and (ws * Kg) A B is small when 6 goes to 
zero, due to the following computation: we have 



{we AB)* Ks{t, x) - {we * Ks) A B{t, x) = j w^{t, X - y)ns{y) A {B{x -y)- B{x)) dy 
hence 

[we AB)* Ks{t, x) - (we * ns) A B{t, x) = j j V ' VS(x - ay) A Wt,e{x - y)Ks{y) dyda. 
To conclude we need to take the Lp' norm in x of that quantity, and Young's inequality yields 

\\{We AB)* Ks{t,-) - {We * Ks) A B{t, ■)\\L2(_n) < Cb\\\ ' I (O) l^e (*) II L2 (H) < C''^II^^O IIl2 

uniformly in time. The result follows for the first equation in (|4.8|) . The second one is of the 
same type, since 

di{we,ib) + d2iwe,2b) = bdiVh'We,h + w^^h ■ "^hb. 
The term {wg^^h " ^hb) * ns is approximated by (w^^h * f^s) ' ^hb exactly as above; to replace 
the term {bdiv^Wg^h) * by ^(div/iti;^ * k^) we write the same type of computation, with 



{bd\VhWe,h) * K5(t, x) - b{d\VhWe,h * H5){t, x) = - J d3We,3it, X - y)Ks{y){b{x -y) - b{x)) dy 
hence, since n_|_ commutes with d^, 

\\{bdiVhWe,h) *KS - b{diVhWe,h * '^5)||L2(R+xn) < CB\\d3We,3\\L^{Il+xn]^ < CbS\\uo\\l2. 

That ends the proof of (|4.8|) . 

Now to end the proof of the proposition, we still need to check (|4.7|) . The idea is to use the 
following estimate, due to the fact that is bounded in Lf{JC) for some compact subspace 
of L^(r2): there is a continuity modulus uj such that 

(4.10) Vy G ri, \\weit,- + y) - Weit,-)\\L2(^Q) < uj{y), uniformly in t and e. 
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Now since ^ = * We^h, we have 

r \ y 

^hwl^h{x) = -^{Vhn){^)we,h{x - y) dy. 
Since Vhi^iy) dy = 0, it follows that 



^hwlhi^) = I ^{\/hK){^){We,h{^ -y)- We,h{x)) dy. 



n 



Then by (lOH) we find that 
where ri{6) goes to zero as 5 goes to zero. The result is proved. 



□ 



Now we are ready to prove the following result. 

Proposition 4.4. Let vP G he a divergence-free vector Geld. For all e > 0, denote 

by Ue a weak solution of i^HJ); and by Ws = Hi_Us its projection onto Ker(L)"'-. Then 

n {we A rot We) in V'{K+ x (O U S)) as e ^ 0. 

Proof. Since the result we are looking for is a weak limit, we can restrict our attention to the 
set O where V5 does not vanish, as in S the result is due to Proposition 13.41 Moreover, to 
prove the result one can restrict our attention to J We A rot We dx3 and it is enough to prove 
that it is proportionnal to V/,6 (up to a small remainder term): taking the scalar product 
with a function in Ker(L) will then yield automatically zero by definition of Ker(L). 

In order to simplify the computations, we shall directly prove the result replacing We^h by ^ 
and Ws^'s by d^W^^. The difference in the two computations is indeed small when 6 is small, 
uniformly in e, exactly as in the proof of l\'6.7\) in the constant case. So writing d^W^^ = 3, 
and dropping the index 5 to simplify, we can perform the following algebraic computations, 
which will prove the result. 

Let us start by recalling that, due to Proposition 14.21 we have 

(4.11) We,h = -^{edtPe,h + re), 

where similarly to Section we have noted 

def vT±TiA ± 
Ps,h = "^h We,3 - W^^^. 

We shall also define 

def „ r^ 

Pe,3 = OlWs,2 - 02Ws,l. 

In ()4.11|) and in the following, the function denotes a remainder term, arbitrarily small in 
the space Lf^ci'^'^ x (O U S)). It follows that 



Wr A rot Wr 



(edtPe,i + re) 


i{edtPe,2 + re) 




22 I. GALLAGHER AND L. SAINT-RAYMOND 

which imphes that 

/ -l{edtPs,2 + re)Ps,3 - d3We,3d3Pe,2 

We A rot We = \{edtPe,l + re)Pe,3 + d^We^sdsPe,! 

\ \{edtPe,2 + re)d3Pe,l - \ {edtPe,l + re)d3Pe,2 

Since the vertical component can be treated exactly as in the constant case, we shall now 
restrict our attention to the first two components. So calling ah the horizontal components 
of that vector field, we have after an integration by parts and using the divergence free 
condition dssWe^s = -<iivhWe,h- 

j ahdx3 = - ^ j {edtPe^h + '^e)Pe,3 dX3 - j pjf^diVhWe,h dX3. 

Now we recall (calling once again generically the small remainder terms) that 

-£dtPe,3 + diVhibWe,h) = re 

so 

diVhWe,h = -^{edtPe,3 - We,h ' ^ hb + re). 

It follows that 

J ahdX3 = j edtPe^hPe,3dX3-^ j p^h^dtPe,3 dX3 + ^ j p^h'We,h-^hbdX3+ J redX3 

^ ~\ j ^^tipthP<^,^)'^^3 + ^ j pi^h'We,h-^hbdx3 + j redx3, 

where now denotes generically the product of by a component of pe- But pe^3 is a 
combination of derivatives of We whereas Pe^h is a combination of components of We- A 
product of the type Pe,hre clearly goes to zero in X''(R+ x (O U S)). For the term rePe,3, one 
uses result (|4.7|) and (|4.9|) to infer that for any subset if of O U 5, 

\\rePe,3\\L^{K) < \\re\\L^{K)\\Pe,3\\L^{K) 

< (csOeil) + C6)\\Whwlh\\L^^K) ^ 0, 

as e followed by 6 go to zero. So from now on will denote generically a term going to zero 
in P'(R+ X (OUcS)). 

Recalling (|4.11j) we get 

j ahdx3 = j edt{Pe^hPe,3) dx3- ^ j Pe^h^dtPe,h ' ^ hb dx3 + j re dx3. 
In particular we have 

j aidX3 = -^ j edt{pe,2Pe,3)dX3-^ j Pe,2£dtPe,\dlbdX3-^ j edt{pl^2d2b)dx3+ j redX3. 
Similarly 

j a2dx3 = ^ j edt{pe,iPe,3)dx3 + ^ J Pe,i£dtPe,2d2bdx3 + -^ J edt{pl^idib) dx3 + j redX3. 
But one can also write 

Ci2dX3 = ^ j £dt{pe,lPe,3+-^pl,ldlb+^Pe,lPe,2d2b)dX3-^ j Pe,2£dtPe,ld2bdX3 + j fgdxs. 
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It follows that up to full derivatives of the type edt, J ah dx^ is equal to 
(4.12) ^--^ j Ps,2£dtPe,idx3^Vhb + j dx^. 

Now the proof is almost finished: we recall that we want to take the projector of the 
term j /\ rot dx^ onto the kernel of L restricted to the set O. Recalling that KerL 

is made of vector fields of the type V^(/3 with V/i6 • = 0, we have obviously 

n ^-^ j Pe,2edtPe,i dx^ Vhb = 0. 



The result is therefore proved for the horizontal component of J We A rot Ws dx^ . The third 
component is identical to the constant case, so the proposition is proved. □ 

This ends the proof of ()4.3() . hence of Lemma l4.ll □ 

In the following we shall denote by {C) the limiting system: 

dtu-i^UAhU + Il{u-Vhu) = in P'(R+ x (O U 5)), U±u = 

u\t=o = nu°. 

The existence of solutions to this system is easy to prove, as it is of the same form as a 
2D Navier-Stokes equation: the only point which might be a problem is that 11 does not in 
general commute with the Friedrich frequency truncation J„, recalled in Section [2| however 
approximating the system by for instance 

will do the job. In any case in the next two sections we shall give precise formulations for the 
solution of (£): in 5 u/j is the unique solution of {NS2D), and outside S it is the solution of 
a heat equation. The uniqueness of u implies in particular that the convergence holds for the 
whole sequence Ue and not only for a subsequence. 

As noted earlier in this section, the third component of Iln for any vector field u is the same 
whether b is constant or not and is simply the vertical average of . It follows that the third 
component of this equation is simply 

dtU3-uAhU3 + Uh-'^hU3 = in R+x(Ou5), U3\t^Q = j U3\-t=Q dx3. 

Now all the work consists in determining Uh. We shall consider separately the vector field 
on O and on 5, which is the object of the two following sections; so in those sections, our 
attention will be restricted to the horizontal component 11^. 



4.2. The averaged equation on O. We shall prove the following result, which yields the 
part of Theorem 121 which lies in O. 

Proposition 4.5. Let u he a vector field satisfying (C) with 

uG L~(R+,H)nL2(R+,iLi) and u £ Ker{L). 
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Then the vector held Uh satishes the foUowing heat equation: 

dtUh - i^nA/iM/j = in R+ x O 

Uh\t=0 = ^Uh,0- 

Proof. The function satisfies 

dtUh - i^UAhUh + n(u/i • VhUh) = 0, 
and since u is in Ker(L), we liave 

- {UAhUh\uh)L2(0) = ~ i^hUh\uh)L2(^o) = \\'^hUh\\l2(^ay 

We note that Uh is equal to zero on the boundary of O, since it is a multiple of V^6. 
So to prove the proposition, the only point we need to check is that 

(4.13) V$ G Ker(L) n H''{0), s > 1, {uh ■ VhUh\^h)L^^o) = 0- 
By definition of Ker(L), we have 

<^h ■ ^hh = Uh ■ ^hb = 0, 

from which we infer that 

(4.14) ^h^Uh = Q. 
Now as in we can write 

n(ll/i • VhUh) = -n(n/j A ToiUh) 

hence 

{uh ■ V/,II/i|^>/i)^2(0) = {uh A ^>/i|rotu/,)^2(o) . 
Finally Identity ()4.14|) yields ()4.13|) . and Proposition 14.51 is proved. □ 

4.3. The 2D Navier-Stokes limit on S. In this section we shall analyse the equation 
satisfied by the limit on 5, that is to say in the regions where 6 is a constant. 

Proposition 4.6. Let u he a vector Geld satisfying (C) with 

uE L~(R+,H)nL2(R+,i7i) and u £ Ker{L). 

Then the vector held Uh satishes the two-dimensional Navier-Stokes equations in S, with 
homogeneous Dirichlet boundary conditions: 

dtUh - ^^hUh + Uh ■ ^hUh = -^hP in R+ x 5 

Uh\t=Q = J Uhfl dX3 

Uh\dS = 0- 

Proof. First let us recall why the equation on is the two-dimensional Navier-Stokes equa- 
tion: we simply consider the weak formulation of the original rotating fluid equations and take 
its limit, by integrating against a test function divergence free and compactly supported 
in R+ X S. The weak formulation is as follows: 

(4.15) / / {-Ue ■ dt^ + uVue ■V<^-Ue^Ue- V$) dxdt = / u°<^\t=o dx. 
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Then taking the hmit as e goes to zero m (|4.15|) yields, due to Lemma 

J (-U • 5j$ + z^Vn • V<I> - n (g) u • V<I>) = J u^^^f^Qclx 

where we have noticed that on 5, we have 

n(n • Vn) = V- (n(g)u). 

Now recalhng that u only depends on the horizontal variable, we deduce the expected equation 
on u, up to the boundary terms. To get the boundary terms, we simply recall that the limit u 
is in Lf^^iYC^ , (Qh)) hence cannot have a jump on the boundary of S. Then we notice 
that dS C do, simply because if x € dS, then V6(x) = 0. So the result follows directly: the 
boundary condition is a homogeneous Dirichlet boundary condition. 

Theorem 131 is proved. 

□ 
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